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Abstract. Global physical properties of random media change qualitatively 
at a percolation threshold, where isolated clusters merge to form one infinite 
connected component. The precise knowledge of percolation thresholds is thus of 
paramount importance. For two dimensional lattice graphs, we use the universal 
scaling form of the cluster size distributions to derive a relation between the 
mean Euler characteristic of the critical percolation patterns and the threshold 
density pc- Prom this relation, we deduce a simple rule to estimate pc, which is 
remarkably accurate. We present some evidence that similar relations might hold 
for continuum percolation and percolation in higher dimensions. 
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Consider a regular d-dimensional lattice where a fraction of sites is selected 
independently with probability p and deemed 'black', with the complementary vertices 
said to be white. The aggregate of these spatial lattice elements forms a random 
pattern, which we may partition into clusters after specifying a neighborhood. This 
simple set-up constitutes the standard model of BernouUian percolation theory, and 
is applied to problems as diverse as transport in disordered media, epidemics, and the 
quark confinement transition in the early universe [U El E] . The central result of this 
theory is the existence of a sharp threshold value < p^ < I in an infinite lattice of 
dimension d > I: When p increases across Pc, a single infinite cluster appears almost 
surely and grows in mass with increasing p beyond Pc [H |B] . 

Research in percolation theory focussed predominantly on the universal critical 
phenomena showing up in the vicinity of the threshold. On the other hand, for 
practical application of percolation concepts, it is the specific and non-universal value 
of Pc which is of primary importance. Exact values of Pc are known only for special 
classes of 2-d lattices |71 [HI IB 1131111 CI]- In all other cases, values of pc are estimated 
numerically with computer simulations, which often are time consuming, in particular 
in three or higher dimensional lattices. 

Here, we investigate the signature of the percolation transition in the Euler 
Characteristic of the spatial pattern formed by the percolating clusters. Its mean 
value per site, provides a topological descriptor, which for a lattice A turns out 

to be an exactly calculable finite polynomial in p. For 2d-lattices the polynomials x{p) 
have one non-trivial zero < po{A) < 1. From the comparison with known threshold 
values, we observe that po(A) gives a tight upper bound for Pc(A) for many lattices, 
but we also find exceptions to this rule-of-thumb. In the case of 3d-lattices, each 
x(p)-polynomial has two distinct nontrivial zeros which are again slightly larger than 
the thresholds values of the two distinct percolation transitions of black and white 
clusters. 

For 2d-lattices, we explain this peculiar ordering of po and Pc using the known 
scaling expression for the critical percolation clusters at Pc- Moreover, this approach 
leads to a surprisingly simple relation which combines via x{p) the specific lattice 
geometry with universal critical percolation features into an accurate parameter-free 
estimate of percolation thresholds of all 2d-lattices considered in this note. Our work 
also applies to bond percolation problems when they are reformulated as the equivalent 
site percolation problem on the covering lattice. 

1. The Euler characteristic in percolation theory 

The percolation transition is a paradigm of a non-thermal phase transition, where the 
local merging of black clusters causes an abrupt change in the large scale-connectivity 
of black vertices. Since the Euler Characteristic (EC) is a prominent descriptor of 
global aspects of spatial patterns, we may expect it to be also a valuable tool in the 
study of the percolation transition. In this section, we introduce the EC descriptively 
and discuss its salient features; a more technical but elementary outline can be found 
in the supplementary notes. 

For the time being we consider planar lattices with cyclic boundary conditions. 
The basic object in site percolation are clusters of vertices, naturally defined by the 
connectivity of the host lattice: Two black vertices belong to the same cluster if 
they are joined by a path of black nearest-neighbors on the lattice. Moreover, each 
configuration of black clusters specifies in a natural way an aggregate of white clusters 
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Figure 1. a) Black clusters partition the white vertices into an aggregate of 
complementary clusters. For the square lattice, white vertices are connected by 
lattice bonds and by diagonal bonds across the faces of the lattice, b) The black 
and white clusters of size 1 and 2 on the square lattice. The white perimeter sites 
of a black cluster are connected (and vice versa). 



with a complementary neighborhood, which is in general distinct from the black one; 
it may be visualized by drawing "matching" bonds diagonally across polygonal faces 
of the host lattice, as illustrated in Figure [TJi for the square lattice. The lattice 
comprising the vertices of the original lattice and all bonds between neighboring white 
vertices is called the matching lattice. The complementary connectivities of black and 
white clusters imply, that the white perimeter vertices of black clusters form closed 
boundaries of holes in white clusters and vice versa, as illustrated for black and white 
clusters of size one and two in Figure [T]3. The percolation thresholds of a lattice A 
and its matching lattice A add up to one [T3] 

P,(A)+p,(A) - 1, (1) 

implying that the black and white clusters are simultaneously critical at p = Pc- 

Let gst denote the number of black clusters per lattice site with a fixed size s and 
a fixed number t of perimeter vertices. The mean density of finite black clusters is 
then given by n{p) — J2s ''^s{p) — J2st P^9stP^ , where p — 1 — p is the density of white 
vertices, comp. Figure [TJd. Correspondingly, the mean density of the complementary 
finite white clusters reads n{p) — X]s"'s(-P) — J2stP'^9stP* ■ their pioneering paper 



[13] on exact percolation thresholds in two dimensions, Sykes and Essam considered 
the difference 

X{p) := n{p) - n{p) = -x{p) (2) 

and found that xip) is a finite polynomial, which they called matching polynomial. 
This observation enabled them to obtain, for instance, the exact critical probability of 
site percolation on the "self-matching" triangular lattice from Xtri{Pc) = at pc = \. 

When specialized to spatial patterns P/y on a planar lattice with N sites, the 
definition of the Euler characteristic, X, reads 

X{Pn) = # clusters of - # holes in Pa,. (3) 

From comparison of this definition with Equation [2j and the observation that 
complementary (finite) white clusters constitute the holes in black clusters (see above), 
we see that the matching polynomial x{p) may be identified with the mean Euler 
characteristic per site (MEC) of the black clusters, x{p) = limAr^oo(;^'^(-P/v))p- 
Equation |2] expresses the fundamental topological invariance of the EC, but the 
representation as the difference of two infinite series is not convenient for practical 
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Figure 2. The MEC of the square lattice (left) and of the triangular lattice 
(right). Note, that po is slightly above pc for the square lattice. The triangular 
lattice is self-matching and Pc = Po = 1/2- 



computation. For that purpose we employ Euler's polyhedral formula, which expresses 
the EC 



(4) 



in terms of the number of black vertices edges joining black vertices #e and 
polygonal faces with black boundary respectively. The mean value x{p) is now 
obtained by a simple local calculation. As an example, consider the square lattice: 
(i) A vertex is black with probability p. (ii) The two vertices bounding an edge are 
black with probability p^, and there are two edges per vertex, (ii) The four vertices 
surrounding a face are black with probability p*. Hence, we find for the square lattice 

Xsgip)^p-2p''+p\ (5) 

The graph of Xsq{p) is shown in Figure [2^. Analogously, one finds for the triangular 
lattice Xtri{p) — p — + 2p'^, which has the above mentioned self-matching property 
x{p) = — x(l ~ P) (see Figure |2|3). Whenever the lattices cells have finite number of 
boundary vertices, x{p) is a finite polynomial. 

The graph of Xsq{p) in Figure [2] is typical for the MEC of 2d-lattices. At small 
values of p, black clusters are finite holes in a single infinite white cluster. As long 
as p is well below pc, the density of holes within the small-sized black clusters is 
negligible, hence xip) is positive and increases with increasing p. On the other hand, 
for p > Pc, and 1 — p <C 1, there is a single infinite black cluster with finite (white) 
holes and thus x{p) < 0, in accordance with Equation [s] In the intermediate range 
of p, which includes Pc, the MEC decreases as the black clusters grow in size and 
merge to generate a single infinite component as p passes the percolation threshold at 
Pc- We see that the typical features of the MEC are governed by the interplay of the 
complementary finite black and white clusters. The percolation transition with the 
singular emergence of an infinite cluster leaves its signature only in the zero crossing 
of x{p) at < Po < 1 with a value of po expected to be comparable with pc- 



2. Percolation thresholds and the zero crossing of the MEC 

The idea that the zero crossing of the MEC should occur near the critical probability 
Pc is plausible and is supported, for instance, by the fact that pq — Pc — 1/2 for site 
percolation on the self-matching triangular lattice, but it calls for a more precise and 
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Figure 3. The percolation threshold pc is slightly below the zero crossing po 
of the MEC whenever pc > 5 ■ This order is reversed, if pc <i 1 ^ apparent in 
the panel on bond percolation. The solution p* of Equation |18| provides a very 
accurate estimate of pc- The deviation \pc — p*| exceeds 0.01 only for very open 
lattices with high percolation thresholds. Lattices are in the order of decreasing 
percolation threshold. For vertex configurations, numerical values of pc and p* 
for the 2-uniform lattices, see Tab. 1 in the supplementary material. Numerical 
estimates of p*''*^ for the Archimedean lattices are taken from | 14| . values for 

pbond j^j,g fj,Qnj [151 [Tel . 

quantitative argument. Here we will first compare po with pc and we shall find that 
Po provides generally a tight upper bound to Pc whenever pc > 1/2. 

2.1. Archimedean lattices 

To begin with, we consider the eleven Archimedean lattices, where all vertices are 
equivalent up to a symmetry operation and all faces are regular polygons. The most 
prominent members of this class are the triagonal, the square, the honeycomb and 
the kagome lattice. Each lattice is uniquely characterized by the number of edges 
Hi of the polygons surrounding a vertex |17l . A lattice of coordination number 
z is therefore conveniently denoted by the symbols (ni , . . . ,nz), where identical 
consecutive polygons are often abbreviated as n°'. For example, a vertex of the 
square lattice is surrounded by four squares. The vertex type of the square lattice is 
therefore (4,4,4,4) or (4'') in the abbreviated notation. Similarly, the a vertex of the 
kagome lattice is surrounded by alternating triangles and hexagons and has vertex 
type (3, 6, 3, 6). With this notation, the MECs of the Archimedean lattices are given 



Topological estimation of percolation thresholds 



6 



by 

i=l ' (1=0 ) 

The percolation thresholds of Archimedean lattices are known to very high precision 
[2]. In Figure |3] we compare Pc to the zero crossing po of For the self-matching 
triangular lattice Po = Pc = 5 ■ Furthermore, a close relation appears to exist between 
Pc and Po even for lattices with pc > 5: Pc is bounded from above by po, with po — Pc 
increasing steadily with pc — 5 ■ Similar observations can be made for the duals of the 
Archimedean lattices, see Fig. 2 of the supplementary material. 

2.2. 2-uniform lattices 

Next, we study the larger class of 2-uniform lattices, which again consist of regular 
polygons but have two distinct vertex types [17 . The two vertex types can 
occur with different abundances, and a 2-uniform lattice is commonly denoted by 
Si{n\, . . . , nl_^) + S2(rii, . . . , ri^^), where Si is the fraction of vertices that are of type 
i. The MECs are given by the straightforward generalization of Equation [6] 

1^=1,2 i=l « p=0 J 

which can obviously be generalized to any finite number of vertex types. We are 
not aware that the percolation thresholds of 2-uniform lattices have been previously 
determined, and we therefore estimated them using an algorithm adopted from [15] . 
Again, po provides a tight upper bound to pc as shown in Figure [3j The vertex 
configurations, as well as po and the estimate of Pc for the 2-uniform lattices are given 
in Tabic 4 of the supplementary material. 

2.3. Bond percolation in two dimensions 

Every bond percolation problem is equivalent to site percolation on the covering 
lattice. The covering lattices, however, are not necessarily planar but decorated 
mosaics. A decorated mosaic is constructed from a planar lattice, where in a subset 
of the faces all diagonal connection have been added (the face is decorated) . A pair of 
lattices where complementary sets of faces are decorated constitute a pair of matching 
lattices in the sense of Equation [l] PJ3j . Calculating the MEC of decorated mosaics is 
slightly more laborious, but a general framework for the calculation has been presented 
in ref. (19J. Using this framework, we calculated the MEC of the covering lattices of 
all Archimedean lattices with vertex (ni, . . . ,nz) 

i 

Comparing numerical estimates of Pc |15l I16j to po confirms p^ > p^ ii Pc > \ 
(Figure|3|, albeit with one notable exception for the lattice 6 with vertex configuration 
(3,4,6,4). For lattices with pc < \ the order of pc and po is reversed, as expected 
from the matching properties of Pc and x(p). 

The relation between pc and po is not restricted to regular lattices but also holds 
for the quasi-periodic Penrose tiling and random tessellations of the plane such as 
Voronoi and Delauny tessellations, see supplementary material for values of Pc and po 
and the polynomials of the MECs. 
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Pdec 

Figure 4. a) A lattice face is decorated, when all diagonal connections across 
the faces have been added to the lattice graph, as illustrated here for three faces 
of the hexagonal lattice. We consider randomly decorated lattices where each 
face is decorated with probability pdec- b) The site percolation threshold Pc{pdec) 
decreases smoothly as the degree of decoration is varied from Pdec = to p^ec = 1- 
The zero crossing po(Pdec) of x{P,Pdec) provides a tight upper to pdPdec) if 
PciPdec) > I S'lid vice versa. The solution of Equation 18 p*{pdec): li^s within 
0.005 oi pc(pdec) with the largest deviations at full or no decoration. 

2.4- Randomly decorated mosaics 

Instead of regular decorated mosaics, we now consider lattices where each face is 
decorated with probability Pdec, as illustrated in Figure |4^. We are not aware 
that this type of percolation process, which bears some similarity to a bond-site 
percolation processes, has been studied before. Our numerical estimates of percolation 
thresholds of a randomly decorated mosaic decreases smoothly from pdPdec = 0) to 
PciPdec = 1), in accord with the containment property |20j . Randomly decorated 
lattices fulfill the statistical matching property Pc{pdec) +Pc(l —Pdec) = 1, from which 
Pc(0.5) = 0.5 follows. The averaging over the diS^erent decoration states of the lattice 
is straightforward and the MEG of randomly decorated mosaics can be calculated in 
the same way as that of regular decorated mosaics. For the hexagonal lattice one finds 

XiPiPdec) = (1 - Pdec)Xhex{p) + (P)- (9) 

From Equation [2] follows the symmetry relation x{p ^ P dec) = ^x{^^P^^^ Pdec)- Hence, 
we have Po{pdec)+Pn{^~Pdec) = 1 in analogy to Pc{.Pdec)+Pc{'^-Pdec) = 1- Our results for 
the randomly decorated hexagonal lattice are shown in Figure [4jD. The zero crossing 
PoiPdec) follows PciPdec) vcry closcly, being a tight upper bound for pdec > h and a 
lower bound otherwise. Similar results can be obtained for other Archimedean lattices 
(data not shown). 

3. The EC of critical percolation and estimation of Pc(A) 

In the previous sections, we saw that Pc(A) - a global property of the lattice - is 
followed rather closely by pq{K), a locally computable quantity. Here, we are going 
to explore the relation of Pc(A) with po(A) in more detail by evoking a generally 
accepted scahng form for the densities Us (p) of large clusters at Pc [H] . These densities 
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determine xip) according to Equation [2] Moreover, they also enter in the sum rules 

13 



(10) 



which express the probability for a particular vertex to be black {p) or white (p) . 

In 2d-lattice graphs only a single critical point exists, so that Poo{Pc) = = 
PooiPc)- At the threshold the scaling ansatz reads 

ns{pc) ^ a{A)s^'^ and ns{pc) ^ a{A)s~'^. (11) 

The non-universal amplitudes a(A) and a(A) account for the particular structure of 
the underlying lattice and its matching partner, whereas the value of the universal 
exponent r = 187/91 is known exactly in two dimensions [4j [22| [23]. 
In order to exploit the scaling hypothesis, we define 



Xsoip) xip) 



So 

E 

s=l 



ns{p) - ns{p)] 



and set 



Likewise, 



XsoiPc) ^{a-a)^s ^. 



s>so 



^S„{P) -^(P-P) 



So 

E 

s=l 



s[nsip) -n^ip)]; 



As„(Pc) 



s>so 



After elimination of the non-universal scaling amplitudes, we arrive at 

C(t, So) 



XsoiPc) 



AsoiPc), 



(12) 
(13) 

(14) 
(15) 

(16) 



C(t~1,so) 

where C(t, sq) = '^2s=ii^ ~^ so)~^ is the Riemann Zeta-function with offset sq. 

The relation ( |16[ ) may be applied as an equality, for instance, (i) to determine a 
value of So from the requirement that the left- and right-hand sides equalize within 
a prescribed accuracy, or (ii) to estimate Pc(A). For the latter purpose, we rewrite 
Equation 16 by substituting the defining expressions (13 151 for Xsoip) and Asg{p). 
The result is 



x{p) = 



C{t, Sq) 

C(t-1,so) 



(2p - 1) - ^ s [ns{p) - ns{p)] 



s=l 



(17) 



So 



+ y[ns{p)-ns{p)]. 



s = l 

The real root, p(so), < p{so) < 1, of this polynomial equation provides an estimate 
for the accuracy of which increases with sq. For the square lattice, the cluster 
numbers gst and gst are known up to s = 12 [21] and we calculated p{so) for 
Sq = 0,...,12. Figure [5] shows how p{so) approaches Pc with increasing sq- The 
Zeta functions ^(t, so) and ^(r — l,so) can be approximated by integrals and their 
ratio evaluates to (r — 2)(r — l)^^(so + 1)^^- In many cases corrections for small 
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Table 1. The deviation of the estimate p{so,A) from Pc(A) decreases rapidly 
when so is increased from zero to two. 



so 


63 


4,82 


4,6, 12 


3,122 


1 /3,4,3,12\ 

2 I 3,122 ) 





-0.01027 


-0.02258 


-0.03184 


-0.04543 


-0.05085 


1 


-0.00553 


-0.00901 


-0.01530 


-0.01240 


-0.02093 


2 


-0.00551 


-0.00283 


-0.00701 


-0.00615 


-0.00937 



clusters are not even necessary, and using sp = ^-^id the approximation for the Zeta 
function yields the simple equation 

X(P) = ^(2P-1), (18) 

with a unique solution p* , < p* < 1 . For all lattices discussed so far, p* is a fairly 
accurate estimate of Pc, as shown in the Figures [3| and Further examples may be 
found in the supplementary note. Among the lattices studied here, the deviation 
of p* from Pc was greatest for the lattices with very high percolation thresholds 
(comp. Figure [3^ and Figure ^jp) which are inhomogeneous on small scales. To test 
whether these deviations are caused by the smallest clusters and holes, we enumerated 
the clusters and holes of size one and two for the hexagonal lattice, the (4, 8^) 
lattice, the (4,6, 12) lattice, the (3, 12^) lattice, and the 2-uniform lattice with vertex 
configuration ^(3, 4, 3, 12) -t- i(3, 12^). The deviation of p(so) from decreases when 
the contributions of the smallest clusters and holes are subtracted, i.e. sq is increased 
from to 2, see Table [ij 

The estimation of the threshold value pc via Equation [l8] will fail, if small clusters 
or small holes are much more abundant at the critical point than an extrapolation of 
the asymptotic law for large clusters would suggest. In particular, this is obvious 
for lattices that contain a substructure which does not contribute to large scale 
connectivity but dominates the density of small clusters. 

§ In many cases, the integral approximation of the Zeta function yields better results than the exact 
ratio, which is probably due to a subtle cancelation of errors. 
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4. Three dimensional lattices 

The combinatorial EC of percolation patterns P on a d-dimensional lattice is given 
by the alternating sum of the numbers of A:-dimensional lattice cells, k — 0,. . . ,d, 
contained in P. Thus, in the case d = 3, the Equation [4] is replaced by 

A'(P) - #« - #e + #/ - #c, (19) 

where #c is the number of black three-dimensional polyhedral lattice cells. The 
resulting MEC of black clusters in the case of site-percolation on a simple cubic (sc) 
lattice, for example, is given by 

X.c(p)=p-3p2 + 3/~p8. (20) 

The black vertices are connected, i.e. they are part of the same cluster, if they are 
joined by a path of black nearest neighbors on the lattice. As in the two dimensional 
case, the black clusters partition the white vertices into an aggregate of clusters with a 
complementary connectivity, such that there is a one-to-one correspondence between 
cavities in black clusters and finite white clusters, as well as between the cavities in 
white clusters and finite black clusters. A white vertex of the simple cubic lattice, 
for example, is connected to all 26 vertices on the boundary of the eight surrounding 
lattice cubes, in contrast to the six neighbors of the black vertices, which correspond 
to lattice bonds. In addition to cavities, black clusters can have handles, i.e. they 
can be homeomorph to a solid torus or objects of higher genus. Each handle of a 
black cluster is pierced by precisely one handle of a white cluster. From these duality 
relations and the 3d analogue of Equation |3] 

X{P) = ^clusters — ^handles -I- ^cavities, (21) 

we see that EC of white clusters is identical to that of black clusters. Hence, MEC of 
white clusters is xiP) — x(l ~ P) with p — I — p. 

The graph of Xsc{p) shown in Figure[6]is typical for the MECs in d = 3, where the 
MECs have two distinct non-trivial zero crossings Pq and 1— p^. In the intermediate 
regime where x{p) < 0; the MEC is dominated by interwoven white and black handles. 
The percolation thresholds of the lattice with black p^ and white connectivity are 
different in general. In the range p'^ < p < 1—p^ a single infinite black cluster coexists 
with a single infinite white cluster. In order to check for a possible link between zero 
crossings of x{p) ^^id thresholds we compare Pq and Pq with simulation values 
and for the sc lattice, face-centered-cubic (fee) lattice and the body-centered-cubic 
(bcc) lattice. The calculation of the MECs for fee and bee lattices is reported in the 
supplementary note. As the Figure [b] already indicates, Pq and Pq are both upper 
bounds for and p^, and they are becoming tighter as the (effective) coordination 
numbers increase. The task to device a threshold estimator based on the above findings 
appears to be more difficult than in the two dimensional case, and it is left for future 
work. 

5. Continuum percolation 

So far, we dealt with the EC of percolating clusters on geometric lattices. Let us 
finally make a few remarks to indicate that the features of the EC induced by the 
percolation thresholds persist in the case of continuum percolation. Consider the 
standard Boolean model where penetrable convex grains are positioned randomly at 
Poisson distributed points in R'^. The grains may be multidispersed having random 
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Figure 6. The graphs of the MEC of the simple cubic (sc) and body-centered 
cubic (bee) lattice. The bcc lattice has the same black and white connectivities 
and is hence invariant to the substitution p = 1 — p. 



Table 2. Numerical values for po E'lid Pc of the cubic lattices for black and 
white connectivities. The bcc lattice has equal black and white connectivities. 
Threshold densities are taken from °) I25I , '') I26| 



lattice 


z 


PO 


Pc 


sc (black) 


6 


0.3940 


0.3116" 


fee (black) 


12 


0.2370 


0.1992" 


bcc 


14 


0.2113 


CITS'" 


fee (white) 


18 


0.1616 


0.136'' 


sc (white) 


26 


0.1139 


0.097'' 



size, shape and orientation. By using results from integral geometry, the MEC of 
patterns formed by clusters of overlapping grains can be calculated exactly [57]. In 
d = 2 the MEC reads 

X2{v) = V{I - V)e^'\ V^po., (22) 
with the grain density p, and the dimensionless ratio / — a and u denote the 
average grain area and perimeter, both assumed to be distributed independently 
from the grain locations. In the case of monodispersed convex grains, / > 1 is an 
isoperimetric ratio. For d~'d one finds 

X3(C)-(l-3/iC+^/2C')e-^ C^pv. (23) 

Here, h = tt, h = 2 ; b, s and v are the averages of the grain mean breadth, 
surface area and volume, respectively. Again, for monodispersed convex grains, Ii > 1, 
/2 > 1. The graphs of x{p) ^^'^ xip)i when plotted as functions of the mean coverage 
< 1 - e"'' < 1 and < 1 - < 1 are similar to the graphs of the corresponding 
lattice MECs: X2{v) has a single zero at rjo = XsiC) two zeros located at 
Co = (1 ± ^1 — ^-§2)- For monodispersed discs 770 = 1, which may be compared 
with the threshold value r/c w 1.12. Monodispersed spheres yield Cq' — 0.377, to be 
compared with the percolation threshold of penetrable spheres, Cc ~ 0.34 [28]. The 
interpretation of the continuum MECs in terms of cluster structures can be carried 
over unchanged from the lattice examples. Thus, for instance, (^q is expected to 
provide a quick estimate for the percolation threshold of the void space defined as 
the set complement of the pattern by the clusters of overlapping grains. 
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6. Summary and Discussion 



Taken together, the equations (17) and (18 1 represent our main result. The relation 



(18 1 connects the universal scaling behaviour of percolating clusters at with their 
specific topological structure, as expressed by the mean Euler characteristic. It 
provides a novel and parameter-free estimate of threshold values for two-dimensional 



lattices, which is remarkably accurate. In the case of self-matching lattices Eq.(18l 
reproduces the exact values Pc — 1/2, and offers an explanation for the numerical 
finding that the zero crossing of the MEC, po, is a tight upper (lower) bound for the 
threshold Pc if Pc> 1/2 {pc < 1/2). 

Over the years, a variety of "universal" approximate formulae for predicting 
percolation thresholds have been devised [HI HH EHl EI] ■ Most of these proposals fit 
an empirical relation with a number of free parameters to a set of known thresholds. 



Recently, Wierman and Naor [32] introduced a list of criteria for the evaluation of such 
formulae. Accordingly, these should (1) be well defined, (2) be easily computable, (3) 
provide values only between and 1, (4) depend only on the adjacency structure of 
the lattice, (5) be accurate, (6) be consistent with the matching relationship, and (7) 
be consistent with the containment principle. 



The Eq.(18l complies with the first six of these requirements, as can be inferred 
from its deduction and from the comparision of our estimates with precise numerical 
threshold values. We have not checked in detail the criterion (7). The claim, that 



our pc-estimation is well defined may perhaps be questioned, since Eq.(18l involves 



an extrapolation of the scaling ansatz to the smallest cluster sizes. However, this 



ad hoc simplification can be removed by going back to the "master equation" (17), 
which provides means to systematically correct for small-scale irregularities, as shown 
in Table (T] 

We mentioned the exact expressions of the mean Euler characteristic for 3d- 
cubic lattices and for the Boolean model of continuum percolation in two and three 
dimensions. The close linking of the zero crossings of these MECs to the respective 
threshold densities appears to persist and underlines once more the role of the Euler 
characteristic as the appropriate concept for describing the topological aspects of 
percolation. But to apply this intriguing fact for the construction of threshold 
estimators for three dimensions and for continuum models remains a challenging 
problem. 
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